Pathwise Solution of a Class of Stochastic Master Equations 
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In this paper we consider an alternative formulation of a class of stochastic wave and master 
equations with scalar noise that are used in quantum optics for modelling open systems and con- 
tinuously monitored systems. The reformulation is obtained by applying J.M.C. Clark's pathwise 
reformulation technique from the theory of classical nonlinear filtering. The pathwise versions of the 
stochastic wave and master equations are defined for all driving paths and depend continuously on 
them. In the case of white noise equations, we derive analogs of Clark's robust approximations. The 
results in this paper may be useful for implementing filters for the continuous monitoring and mea- 
surement feedback control of quantum systems, and for developing new types of numerical methods 
for unravelling master equations. The main ideas are illustrated by an example. 
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I. INTRODUCTION 

In quantum optics stochastic wave and master equa- 
tions arise in the study of open systems and continu- 
ous measurement, see, e.g. 0, 0, H, 0, and the 
many references cited in these works. These equations 
are driven by stochastic inputs, typically white noise 
(Wiener process), representing photocurrent, or Poisson 
jumps, representing photon counts, and involve stochas- 
tic integrals — they are stochastic differential equations 
(SDEs), 0,0- These SDEs can be solved or approxi- 
mated numerically either for use in simulating open sys- 
tem dynamics, or for updating conditional states (the 
topic of numerical approximation of SDEs is well docu- 
mented ^ Chapter 10]). However, it is important to 
keep in mind that these are idealized models (e.g. the 
Wiener process is highly irregular, in fact nowhere dif- 
ferentiable with probability one) , and the models may be 
used in conjunction with real data. Hence it is of interest 
to consider the robustness of stochastic wave and master 
equations from a practical point of view. 

Statistical robustness was considered in 1978 by 
J.M.C. Clark in the context of classical nonlinear fil- 
tering. The theory of nonlinear filtering is an important 
and well documented part of the systems and control, 
communications, signal processing and probabi lity and 
statistics literature. It is well known (see, e.g. [lOl lU\) 
that the solution to the nonlinear filtering problem, say 
for a diffusion state process observed in white noise, is 
given in terms of the conditional distribution which solves 
a measure- valued stochastic differential equation (analo- 
gous to the stochastic mater equation) . The correspond- 
ing equations for the conditional density is a stochas- 
tic partial differential equation. Clark drew attention 
to the disadvantages of stochastic integral representa- 
tions of nonlinear filters from a practical point of view. 
These disadvantages concerned lack of statistical robust- 



ness and the inability to cope with the range of measure- 
ment data (driving process) that can arise in practice. 
Clark addressed these issues by providing a reformula- 
tion of the nonlinear filtering equations that does not 
involve stochastic integrals. Clark's so-called pathwise 
solution defines a version of the conditional distribution 
(or density) that is defined for all possible measurement 
data and is a continuous function of the measurement 
data, thereby providing important robustness qualities. 
Clark also provided numerical approximations to the re- 
formulated nonlinear filters which inherit the robustness 
characteristics. For further details, see 0,^3' ^'^d re- 
lated matters [T3| . 

In this paper we give reformulations of the stochas- 
tic wave and master equations with scalar noise that do 
not involve stochastic integrals, analogous to the clas- 
sical pathwise versions of nonlinear filters proposed by 
Clark. The reformulated equations are ordinary differ- 
ential equations where the driving process enters as a 
parameter — not via a stochastic integral. This reformu- 
lation may be useful for implementing filters for the con- 
tinuous monitoring and measurement feedback control 
of quantum systems, and for developing new types of nu- 
merical methods for "unravelling" master equations. 

The paper is organized as follows. In Section ^1 we 
describe the stochastic wave and master equations to be 
considered, and provide some motivation and background 
information. Then in Section IlIII the pathwise solution 
and robust approximation for quantum diffusion case are 
presented, and we give an example to illustrate the solu- 
tion in the context of an imperfectly observed two-level 
atom continuously monitored by homodyne photodetec- 
tion. Section IIVI contains the formulation for quantum 
jump case with some brief comments. Some calculations 
and the proof of a continuity result are provided in the 
Appendices. Some of the results in this paper were an- 
nounced in the conference paper |14| . 
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II. STOCHASTIC WAVE AND MASTER 
EQUATIONS 

A. Background 

We recall (see, e.g. [H Chapter 2], [3) that an 
isolated quantum system is described by a (pure) state 
£ H (Dirac bra-ket notation), where H is a com- 
plex Hilbert space, with time evolution governed by the 
Schrodinger equation 

th-\,Pt) = H\,P,) (1) 

where h = h/2Tr and h is Planck's constant, and _ff is a 
Hamiltonian operator. In what follows we use units such 
that h= 1. 

However, when a quantum system is interacting with 
an external environment, the interactions must be taken 
into account. In the open systems literature (see, e.g. 0, 
Chapter 5.4], Chapter 6]), the Schrodinger equation 
(Q) is replaced by a master equation, which takes the 
form 

p = -i[H,p]+V[L]p; (2) 

here p is the density operator and the superoperator V 
is defined for any operator c by 

The system operator L is used in the modelling of the 
interaction. Note that when L = and p = \tp){il'\, the 
master equation (j^J reduces to the Schrodinger equation 

m- 

A common method for solving the master equation ^ 
(via "unravelling" ) is to first solve a stochastic wave equa- 
tion and then to average. For instance, one could solve 
the linear equation 

d\4!) + K\iP)dt ^ L\iP)dy (3) 

for an unnormalized state \ip), where 

K = iH +^L^L, (4) 

and y{t) is a Wiener process, and then average 

pit)^n\m)m)\i 

where 

Vwwi*(*)> 

Such a procedure is computationally advantageous since 
the wave function contains fewer components than the 
density operator. 



Stochastic wave equations of the form (jJJ and related 
stochastic master equations also arise when quantum sys- 
tems are continuously monitored P, Chapter 11], 0. To 
motivate this, we recall that an ideal measurement of the 
system is characterized by a self-adjoint operator A on H. 
In the simple case that A has a discrete non-degenerate 
spectrum {ui} C R, the possible outcomes of a measure- 
ment are the eigenvalues a^. The outcome is random, 
where Ui occurs with probability 

when in state (assumed normalized: \tp\'^ = 

tr[\^p){ip\] = 1 ). Here, \ai) denotes the orthonormal 
eigenvector of A corresponding to the eigenvalue a^. Af- 
ter the measurement, there is a collapse of the state to a 
new state 

|V:) = |a.)(a,|^)/V^. (5) 

The state is the conditional state given the mea- 
surement outcome a^. Consequently, when a quantum 
system is measured, the deterministic evolution given by 
the Schrodinger equation ^ must be augmented by a 
stochastic transition, e.g. ((SJ). The continuous measure- 
ment of quantum systems can be regarded in terms of 
a sequence of measurements of infinitesimal strength (in 
contrast to the possibly large jump in which accumu- 
late in the limit to provide the conditional information 
and the evolution of conditional states as in Q and re- 
lated stochastic master equations (see below), P, 
Chapter 11]. 

In this paper we consider two kinds of stochastic mas- 
ter equation corresponding to the two standard types of 
stochastic integrator with independent increments: the 
standard Brownian motion (diffusion case) and the stan- 
dard Poisson type (jump case). 

B. Quantum Diffusion 

We consider a stochastic master equation (SME) 
dpt = [LptL^ ~ Kpt - ptK^]dt 

+ - [Lpt + ptL^ ~ ptMp^] dvt. (6) 

In ©, the operator K is defined by Q. In the case of 
continuous measurements, the parameter k > 1 is related 
to a measurement efficiency parameter < 77 < 1 via k = 
1/y^ corresponding to imperfect or noisy measurement; 
here perfect measurement corresponds to = 1. 

The SME jnj is driven by real valued white noise z>t, 
represented in ^ by an Ito-sense stochastic integral with 
respect to a standard Brownian motion (Wiener process) 
ut, sometimes called an innovations process. This process 
is related to a real valued process j/t, which we call the 
measurement process, by 

dyt = Mp^dt + ndvt, (7) 
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where Mp={L + Lt)^ = tr{(L + L'^)p}. 

If pt denotes the expected value of pt, then pt solves 
the master equation |2Il. Note that in the case L — Q 
(no measurement or interaction) and initial pure state 
Po = \i'o){'4>Q\, the state is pure for all t, pt = |V't)('0t|, 
and © reduces to the Schrodinger equation for jV-'t) 

The SME © is nonlinear in pt (due to the term 
ptMp^), and the solution pt of © is normalized for all t: 
tr Pt — 1. We find it convenient to work with an unnor- 
malized version pt, defined by 



where 



At 



exp 



Pt = Atpt 



^ Jo 



(8) 



It can be checked using Ito's rule (see, e.g. 0, Chapters 
12 and 18], TT', Chapters 6 and 7]) that pt solves the 
following linear stochastic equation: 



dpi 



= [LptL^-Kpt-ptK^]dt+^ [Lpt + ptL^ dyt. (9) 



Note that this unnormalized SME is simpler and in bilin- 
ear stochastic form driven by the measurement process 
dyt- The unnormalized density operator pt can be nor- 
malized by simply dividing by its trace (Aj — tr(pt)). 
Equation (O is analogous to the Duncan-Mortensen- 
Zakai equation of nonlinear filtering f [loL Chapter 18], 
[H Chapter 7]), and is also known in the quantum 
physics literature (,19j, Section 4], p^. Section 4.2.2]). 

Note that in the case k — 1 and initial pure states, 
the unnormalized SME © reduces to the unnormal- 
ized stochastic Schrodinger equation |(2Jl; see |22t Section 
4.2.1], ill Chapter 11]. 



counting observation process dNt which is a real random 
variable satisfying 



EldNf 



dN^ = dNt 



r]Xtj:{Jpt)dt 



(12) 



The observation process dNt only gives the value of either 
zero or one (corresponding to the counting increment), 
and is the representation of the standard Poisson process 
with intensity 7]\tv{J^pt) dt. Here we assume that the 
observation process has efficiency < 77 < 1 as the rep- 
resentation of imperfect or erroneous counting process. 

The SME together with lO simply tell us the 
behavior of quantum jump that in the increment of time 
dt the system jumps via superoperator J with prob- 
ability Pj = E[dNt] or smoothly evolves via the first 
bracket term of RHS of ifTUI) with probability Ps = l — Pj- 
Note that the average or expectation value of pt in Hl()|l 
obeys master equation |2Jl by considering the operator 
L = Ai/2(c-/). 

One can check for perfect counting process (ry = 1) 
and initial pure states that the jump SME H1U|I reduces 
to the normalized jump stochastic Schrodinger equation 
(see e.g. H 113, EH) 



c 



-I 



iE 
\^t)dNt 



IV-t) dt 



(13) 



The jump SME Hl()|l and 113|l are normalized but non- 
linear so again we work with unnormalized version pt, 
defined by 

Pt = Atpt 



C. Quantum Jumps 

Another type of continuous measurement is described 
by the counting observation which give rise to a jump 
stochastic master equation 

dpt = [-Gpt- ptG^ + {l-ri)XJpt+ifi\ptiy:{Jpt)\dt 
Jpt 



where 



tr(Jpt) 



Pt 



dNt 



The operator G is defined by 



A 



G = -C'^C + iE 



(10) 



(11) 



where C is some Schrodinger picture system operator, E 
is energy operator related to system Hamiltonian hy H ~ 
E + ^{C - C''), \ > is related to the intensity of the 
standard Poisson process and J'pt = CptC^ is defined as 
the jump superoperator. The SME H1U|) is driven by the 



At = 1 + / AMJps - Ps) [dNs - r]X ds] . (14) 
Jo 



And by Ito's rule for jump process, the unnormalized pt 
solves the following linear stochastic jump equation 



dpt = [-Gpt - ptG^ + {1- r])\Jpt + r]Xpt]dt 
+ [Jpt-pt]dNt. 



(15) 



We see that the unnormalized jump SME H15|) re- 
duces to an unnormalized and linear jump stochastic 
Schrodinger equation 

dli^t) = [Hl- C^C) - iE] \iit) dt+[C- I] \^t) dNt 

for the case of perfect measurement and pure initial 
states. 
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III. DIFFUSION CASE 



A. Pathwise Solution 



We follow Clark's approach Q to obtain a pathwise 
solution to the SME ©. Let 



L l? 



(17) 



Let pt be a solution of the unnormalized stochastic mas- 
ter equation and define an unnormalized state rt by 



rt = AtptA. 



(18) 



Then, as shown in Appendix ^ rt solves the pathwise 
master equation 

rt = LrtL^ [l - I/kI'] - AKA^'rt - rt{Al)-' Aj. 

(19) 

Conversely, solutions pt, pt to the stochastic master equa- 
tions ©, ® can be obtained from a solution rt to l(TI?|) 
via the formulas 



Pt ^ A-'rtiAl)-\ Pt 



Aj'rtjAl)-^ 
tr[4-V,(4)-i] 



(20) 



This result is analogous to the classical result for non- 
linear filterine; % Theorems 4 and 6]. It is important 
to note that the pathwise equation 1)19(1 does not involve 
stochastic integrals; the measurement path enters as a 
parameter in an ordinary equation. In particular, the 
versions of the solutions to the stochastic master equa- 
tions ijnji, ® defined by ifT^ . if^ are defined for all con- 
tinuous observations paths, not just for a set of paths of 
full Wiener measure as is the case for solutions obtained 
directly due to the stochastic integrals in JBJ, ©. 

We next make explicit the continuous dependence on 
the measurement paths. We make use of the supremum 
norm 

II / IIt- sup \ft\ 

0<t<T 

for a continuous vector or matrix valued function, and | • | 
denotes the appropriate Euclidean norm. 

Let H be finite dimensional. Then, as shown in Ap- 
pendix^ solutions Pt, Pt to the stochastic master equa- 
tions Q , © defined by (|20|l are locally Lipschitz continu- 
ous functions of the observation trajectories. This means 
that if ul and Ut are two continuous observation trajec- 
tories on [0, r], with corresponding solutions p}, p^and 
"pI, p1, respectively, then there exists a positive constant 
C = C{\\ ||t,|! ||t,T) such that 

IIp'-p'IIt < C\\y'-y^\\T, and 
Wp'-p'Wt < C\\y'-y'\\T. (21) 



We note that in the case of perfect measurement and 
initial pure states, the pathwise SME lfTO|) reduces to a 
pathwise Schrodinger equation: 



h) =-AtKAi^\cf>t 



(22) 



This equation is also defined for all continuous observa- 
tion paths, and depends continuously on them. 



B. Robust Approximation 

In this section we use the pathwise master equation 
(|19|l to derive an approximation to the stochastic master 
equation (|SJ). We will restrict our attention to the case 
of a finite dimensional underlying Hilbert space, and we 
employ a simple implicit Euler scheme to illustrate the 
ideas. In 0, the corresponding approximations for non- 
linear filters were called robust approximations. 

Fix the interval between sampling times A—tn — ^n-i. 
A reasonable implicit Euler approximation for (|19() is 



rtr 



rtti+Lrt^L^ [l - l/«2] A 



i„ A , „ A(^t yiK^l ]A. 



- [Au^KAtWtn +rt 
Multiplying both sides by A^J^{. . .}{AI^)~'^ gives 

Pt^ = Ar>t„-iPtli4„_,(4j-' 

+ Lpt^L^ [1 - l/«2] A - [Kpt^ + pt^K^]A. 
Defining Ay„ = yt„ - Vt„^i, we write 

Au^t^-i = exp <^ — Ay„ - — A \ = £{Ayn) 

<_^(4j-i = e^J^^Ay,^^^^A^E{AyJ. 

We obtain the implicit robust approximation for the un- 
normalized SME expressed as a matrix equation 

Apt^ + pt^B - Cpt^V = £iAy^)ptti£i^yny, (23) 
where 

A = [I + KA] 
B = K^A 
C = L 

V = [1-1/^2] A. 

The matrix equation (|23|) can be solved explicitly by 
rearranging elements of the matrix. Suppose A = 
[^i|^2| • ■ • l^n] is an nxn matrix, A\ . . . A^ are col- 
umn vectors of A. We define operator Vec(^) = 
[Ai|^2| ■ • • l^n]"^- Thus, Vec operator transforms nxn 
matrix to (nn) x 1 matrix. 

Suppose B and X are also nxn matrices, then 



N<ic{AXB\ = {B^ ® y^lVecA", 



(24) 
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where (E) denotes Kronecker product. If B is a complex 
matrix, then B"'" is simply transposing B without conju- 
gating it. 

Transform H23|l by Vec operator 

VeciApt^) + Vecipt^B) - Vec{Cpt^V) = (25) 
Vec(£:(Ay„)Ptlif(Aj/„)^) 
[(/ ®A) + [B^ {V^ ® C)] Vec(pt^) = (26) 

Vec{£{Ay.^) Ptt iSi^Vn)^), 

yields 

Vec(p,^^)= [iI®A) + {B^®I)-iV^®C)]'' ,27^ 



Writing H27|l in the symbolic form 

Ptn = r(A?/„)pt^_i, 



(28) 



we obtain an approximation to the solution of the un- 
normalized stochastic master equation 0. By normal- 
ization we have 



Ptr 



T{Ay„)ptti 
tr[r(Ay„)p,^_i 



(29) 



an approximation to the stochastic master equation jnj. 

Note that F in H28|l can be seen as common recursive 
filtering solution incorporating two steps, prediction and 
update or correction. The prediction step utilizes knowl- 
edge in yt histories {yto„-i) via Ptn-i ^-^^ then the re- 
sult is updated by the current measurement information 
available in Ay„. 

A significant result concerning robust approximations 
0, Theorem 7] is that the convergence of the approxi- 
mation pt^ to the exact solution pt^ is pathwise for all 
observation trajectories. Indeed, the following inequality 
can be proven as in Theorem 7]: there exists a con- 
tinuous function fc(-) such that for all A > and n with 

0<tn<T, 

\pt^{y) - PtAy)\ < H\\ v \\t){a + wy{A)) (30) 

where 

?i;3^(A) = max{|j/(si) - ?/(s2)| : 0<si,S2<A}. 

This should be compared with other discrete approxi- 
mations to solutions for SDEs, 0, 0, Chapter 10]; for 
example, as discussed in "P, Section 4] , a direct Euler ap- 
proximation of converges "almost surely" , but for dif- 
ferentiable observation trajectories y, the Euler approx- 
imation converges to a limit which is not the same as 
pt{y) given by (j^OJ). 



is C^, the two-dimensional complex vector space, whose 
elements are called qhits in quantum computing. Let |0) 
and |1) denote basis vectors corresponding to ground and 
excited states, respectively. We use the following Pauli 
matrices to represent operators for this system: 



a. = |0)(l| + |l)(0| = 
a,=z|0>(l|-z|l>(0| 



Il)(l| 



|0)(0| 



1 

1 

-i 

1 

1 0' 

-1 



= |0)(1| = -;^{(^x ~ i(ry) = 





1 



(31) 



Here cr is a system (lowering) operator. Any state p on 
can be represented in terms of the Bloch vector (x,y, z) 
([H Chapter 2]): 



ycTy 



1+z x—iy 
x+iy 1—z 



(32) 



where x^ + 
is given by 



< 1. The Hamiltonian of the system 



a A 1 

H — —a-r H (T^ = — 

2 2 2 



A 

a 



a 

-A 



(33) 



where a is the Rabi frequency, A is the atomic frequency 
minus the classical field frequency. 

The two level atom is coupled to an optical field which 
is continuously monitored by homodyne detection (see, 
e.g. 0, Section 8.7], 0, Section II.C]). The output of 
the detector is a current /c(i) whose mean value is pro- 
portional to the expected field quadrature tr[X^p] deter- 
mined by the phase angle ip of the local oscillator. Here, 



1 



lip 



For if = we are interested in measuring the x- 
quadrature 2Xq = ax, while for (p — 7r/2 we are in- 
terested in measuring the ^/-quadrature 2X^^/2 ~ ~'^y 
Variations about the mean are called quantum noise, a 
key feature of quantum optical systems. 

The stochastic master equation for this setup is equa- 
tion © with H given by H33(l and 



(34) 



C. Example 

In this section, we apply the discrete approximation 
derived in Subsection IIII Bl to a two level atom continu- 
ously monitored by homodyne photodetection 2^ Sec- 
tion HFC]. In this example, the underlying Hilbert space 



where 7 is the spontaneous emission rate. In this case 
K = ■^o'^o' + iH . The corresponding measurement equa- 
tion is {Tj), where Ic{t) — yt is the homodyne photocur- 
rent. The quantum noise {l/y^)i>t is white with variance 

The approximation (|29|) was implemented for this ex- 
ample with an assumed detection efficiency of 77 = 85% 
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(a low value), to compare with the results of [22, sec- 
tion III.C] which used different methods and considered 
the case of perfect measurement efficiency rj — 1. The 
simulation was carried out as follows: 

• Set 7=1. All other parameters are based on 7 unit, 
A=0, Q^="^7- Simulations are conducted with two 
values of ip=0 and ip=7r/2. 

• Time step A=0.01/7, time length T^25/j corre- 
spond to simulation length n=2500. 

• Simulations arc done for single ensemble and 
A^=1000 ensembles. 



• Set pure state initial condition of |V-'o 

1 1 



= M> such 



thatpo = |i/'o)(V'o| 



. 1 



1 1 



This corresponds to 



initial Bloch vector {x,y,z)o — (1,0,0). 

• We compute recursively pt via (|29|l . The new mea- 
surement data generated by 

Ay„ = tr{(i + Lt)pt^_ JA + kAi^„, 

where Aj/„ is an independent identically dis- 
tributed Gaussian sequence with mean zero and 
variance A. 

• We obtain the corresponding Bloch vector 

{Xtn,ytn, Ztn) by USiug (|22J). 

In spite of the poor measurement efficiency, one can 
still infer important physical information about the sys- 
tem as described in [24 section III.C]. Indeed, in terms of 
the Bloch vector {xt,yt, zt), the homodync photocurrcnt 



is 



Ic{t) = y/j[xt COS tp-yt sin ip] + 



KVt. 



When the local oscillator is in phase with the driving field 
((/3=0), the deterministic part of measurement is propor- 
tional to {(Jx)- This measurement seems to drive the sys- 
tem into an eigenstate of ax- This is shown in FigQ] Al- 
ternatively, one can see from the steady state ensembles 
simulation Fig|Sl that the atom states are concentrated 
near a:; = -1-1 and x = —1. 

In contrast, measuring the quadrature with </? ~ 7r/2 
will eventually force the atom states into the eigenstates 
of Uy. The states are spinning around the sphere toward 
z ~ +1 and z = —\ due to the driving Hamiltonian, 
FigEl FigH 

The effect of imperfect measurement can be seen 
clearly from these results. The Bloch vectors are not 
confined to the surface of the unit sphere, thus the sys- 
tem is in mixed states. This shows that the imperfect 
measurements can cause loss of information about the 
quantum system, but nevertheless the information is con- 
sistent with the perfect case ^2 section III.C]. 




500 2000 
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1500 2000 



FIG. 1: The evolution of one ensemble Bloch vector with <^=0 




FIG. 2: The evolution of one ensemble Bloch vector with 
V3=7r/2 



IV. JUMP CASE 

To derive a pathwise solution for the jump SME I^UJl, 
we choose (see for example (2^) 

^ = C-^', (35) 

and again define an unnormalized state rt by 

n = AtptAl (36) 

where pt is a solution of the unnormalized jump stochas- 
tic master equation (|15|) . Then rt solves the pathwise 
master equation 

n = -AtGAi\t - rt{A\)-^G^A\ -f (1 - ri)XJrt + TyAr*. 

(37) 

Moreover, solutions pt, pt to the stochastic master equa- 
tions lfTU|l . l(T^ can be obtained from a solution rt to l|S7|l 
via the formulas 



pt = Ai\t{A\r\ Pt 



1^.^ z|t\-l 



At'rtjA]) 
tv[A-;\t{A\)-^ 



(38) 



In the case of perfect measurement and initial pure 
states, the jump SME H37|) admits the jump pathwise 
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The Bloch Sphere 
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and enjoy continuity properties. These robustness char- 
acteristic would be useful when applied to quantum fil- 
tering problem such as in quantum feedback control. 

The results we have established are valid for scalar 
measurements, but can easily be generalized to the case 
of multiple measurements provided the interaction oper- 
ators commute. It is not known if the results generalize 
in case the interaction operators do not commute. 
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FIG. 3: The steady state point of 1000 ensembles Bloch vector 
on the Bloch sphere with <^=0 



The Bloch Sphere 
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FIG. 4: The steady state point of 1000 ensembles Bloch vector 
on the Bloch sphere with </p=7r/2 

Schrodinger equation 

100 = I^^AtGA^' + ^I^ (39) 

The pathwise solutions H37|l and H39|l appear as ordi- 
nary equation without stochastic integrals in terms of 
dNt; the counting observation result Nt enters as a pa- 
rameter in At- These solutions are also defined for all 
continuous counting observation paths. 

However, we note that for the case of the jump SME, 
the pathwise solution might not be so useful for computa- 
tion, and existing techniques (see, e.g. P, Chapters 11]) 
may be preferable. 

V. CONCLUSION 

In this paper we have proposed the pathwise reformu- 
lation of the stochastic master equations for two cases; 
quantum diffusion and jump. These reformulations pro- 
vide solution that is defined for all measurement paths 



APPENDIX A: PROOF OF PATHWISE 
EQUATIONS 

In this appendix we prove the assertions of sections 
nil Al and IIVI concerning the pathwise equations. 

We first consider the diffusion case and verify H19|) and 
(|20|l . The calculations are simple but needs frequent and 
careful use of Ito's rule. Differentiate (|18|l . 

dn = dAt.ptAl + Atdpt.Al + AtPfdA\ + AtdpfdA\ 
+ dAt.pt.dAl + dAt.dpt.Aj + dAt.dpt.dAl(Al) 

and lfT7|l . 

dAt , dAt , 1 d^At , 2 
dAt = -d^dy, + —dt+-^dy, 

= + 
L 

= -At^dyt + At^dt. (A2) 
Similarly, 

dAl = -Al^dyt + Aj^^dt. (A3) 

Put 10, HA2p . (|A3p together into (jXTJ, and carefully us- 
ing Ito's rule, we obtain H19|l . 

In the jump case we prove (pTzIl and analogously 
to the diffusion case. We need to use Ito's rule for dNt, 
i.e. dN^ ~ dNt and dNtdt — 0. Any higher order of 
differential involving dt is zero. Differentiate (|36|) . 

dn = dAt.ptA\ + Atdpt.Al + Atpt-dA\ + Atdpt.dA\ 
+dAt.pfdA\ + dAt.dpt.A\ + dAt.dpt.dA\,{A^) 

and lfT7|> . 

dAt = At (£(-1)"^) dNt 

= At{C-^ -I)dNt. (A5) 
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The infinite series arise due to the Ito's rule that higher 
orders of dNf do not vanish. Similarly, 

dA\ ^ A\{{C^)-^ - I) dNt (A6) 

Put (dSJ, (|X5|) . (|X6|) together into (|X4|) . and carefully 
using Ito's rule, we obtain (|37|l . 

APPENDIX B: CONTINUITY PROOF 



Here Ci and C2 depend on |1 |1t, |1 J/^ ||t and T. By 
Gronwall's lemma, we have 

\rl-rj\ < Ci||yi-2/2|iTexp|^*C2ds| 

< C \\y'-r/\\T 
for all < t < r, which implies I|B1() as required. 



We provide the proof of the continuity result (|21|l of 
section lill Al 

In view of H20() . it is enough to verify that 



r 



2 



\t< C \\ - \\t ■ (Bl) 



Let Equation (fH^ defined for all t e [0, T], and we 
rewrite it in term of matrices TZt — Ti-iyt) and A4 such 
that 

ft = ntrt + rtUl +Mirt). (B2) 

Given observation records y} , y^ we have the correspond- 
ing pathwise solutions — rt{y}), = ftiUt) and ma- 
trices n\=n{y]), 7^2 = 7^(y2). 
It follows 



Jo 

+ MirD- Mirl)\ds 

+M{rl)-My)\ds 

< f\\TZl-TZl\\rl\ + \rl\\nl^-nf\}ds 
Jo 

+ tm\\rl-r'A + \rl-r'J\nm}ds 
Jo 

+ I {\M{rl) ~ M{rl)\] ds 
Jo 

< C,\\y^-y^\\T+ fc^\r]-rl\ds 

Jo 
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